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Observations of the Global Positioning System (GPS) will enable a reduced- 
dynamic technique for achieving subdecimeter orbit determination of Earth-orbiting 
satellites. With this technique , information on the transition between satellite states 
at different observing times is furnished by both a formal dynamic model and ob- 
served satellite positional change (which is inferred kinematically from continuous 
GPS carrier-phase data). The relative weighting of dynamic and kinematic in- 
formation can be freely varied. Covariance studies show that in situations where 
observing geometry is poor and the dynamic model is good , the model dominates 
determination of the state transition ; where the dynamic model is poor and the 
geometry strong, carrier phase governs the determination of the transition. When 
neither kinematic nor dynamic information is clearly superior, the reduced-dynamic 
combination of the two can substantially improve the orbit-determination solution. 
Guidelines are given here for selecting a near-optimal weighting for the reduced- 
dynamic solution, and sensitivity of solution accuracy to this weighting is examined. 


I. Introduction 

The Defense Department's Global Positioning System 
(GPS) will be in full operation, with at least 18 satel- 
lites in orbit, by the end of 1992. The two precise radio- 
metric data types available from GPS, P-code pseudo- 
range and continuous carrier phase, will provide accurate 
positioning for users on the Earth’s surface and in Earth 
orbit [1]. Advanced differential techniques incorporating 
GPS data from a global network of ground reference sites 
now promise to provide the subdecimeter orbit determi- 
nation accuracy being sought for a growing number of 
scientific remote-sensing satellites [2-4]. One such satel- 
lite is the U.S. -France Ocean Topography Experiment, 
TOP EX/POSEIDON, to be launched in June of 1992. 
TOPEX, which will fly in a circular orbit at an altitude of 
1,336 km [5], has a formal accuracy requirement of 13 cm 


for the continuous determination of its geocentric altitude; 
in fact, TOPEX ocean science would benefit from an al- 
titude accuracy comparable to the 2.5-cm precision of its 
radar altimeter. A number of proposed international mis- 
sions for the mid- to late 1990s with high-Earth elliptically 
orbiting radio interferometric platforms would also bene- 
fit from a decimeter-level orbit determination capability. 
The Deep Space Network has responsibility for GPS-based 
tracking of TOPEX and for tracking of other satellites in 
high-Earth orbit. 

Differential GPS tracking can take many forms. The 
simplest is relative point positioning using instantaneous 
differential pseudorange measurements to four or more 
GPS satellites. The accuracy of this geometric relative 
positioning is limited primarily by measurement noise and 


13 


GPS ephemeris error, magnified by position dilution of 
precision (PDOP), which is related to observing geometry 
[1]. For a low-Earth orbiter, instantaneous differential po- 
sitioning accuracy is expected to approach one meter [6]. 
In applications where more accurate orbit knowledge is 
needed, a long-arc dynamic solution may be suitable. With 
this approach, the satellite dynamics are constrained by 
physical models, and noisy instantaneous measurements 
obtained over a period of time are combined to improve 
precision and to yield greater information on the user state 
at a single epoch [2,3]. For greatest accuracy, all GPS or- 
bits and some ground receiver positions are also adjusted; 
all solutions are obtained in a reference frame established 
by a small set of fixed (unadjusted) ground sites. In the 
dynamic solution, the transition from satellite states at 
different measurement times to the state at the solution 
epoch is furnished by integration of the equations of mo- 
tion, which are governed by the forces (dynamics) acting 
on the satellites over the time of interest. Any mismod- 
eling of these dynamics will result in systematic errors in 
the state solution — errors which tend to grow as the data 
arc length increases. 

In [4], a long-arc non-dynamic (or kinematic) tracking 
technique is proposed. With this technique, the instan- 
taneous user satellite positions are again determined by 
differential GPS pseudorange measurements; however, in- 
formation about the transition between positions at differ- 
ent times is furnished by the satellite positional change as 
inferred from observations of continuous differential GPS 
carrier-phase data. These data enable many point position 
solutions to be smoothed over long data arcs. Since the 
GPS carrier phase can be measured extremely accurately 
(to 1 mm or better in I sec), with favorable observing 
geometry it can provide nearly ideal state-transition infor- 
mation. And because the kinematic solution is fundamen- 
tally geometric, the solution is free of dynamic modeling 
errors. In exchange, however, the solution has a high sen- 
sitivity to the continuously changing observing geometry. 
To maintain decimeter orbit accuracy, strong observing 
geometry must be ensured by providing sufficient receiv- 
ing channels on the orbiter and a sound global network of 
ground receivers [4]. 

Kinematic tracking discards dynamic orbit models and 
the associated information entirely. That, indeed, is one of 
its principal attractions. Not only is sensitivity to model 
errors eliminated, but the complexity of the solution pro- 
cess is greatly reduced. Nevertheless, one can expect in 
general that an optimal synthesis of dynamic and kine- 
matic information will offer advantages. In the reduced- 
dynamic technique, first proposed in [7], both the dynamic 
and kinematic methods of state-transition determination 


are used with carefully selected relative weighting. The 
weight on the dynamic information is controlled through 
adjusting a set of three process-noise parameters repre- 
senting a fictitious 3-D force on the user satellite. These 
parameters include the a priori uncertainty cr 0 , the steady- 
state uncertainty <j , and the correlation time r. 

Appending process-noise parameters to the satellite 
force model is attractive since, although the fictitious force 
is piecewise constant and therefore discontinuous between 
batches, the satellite state components, derived by inte- 
grating the noise parameters with the force model, remain 
continuous. This fits very naturally with the traditional 
dynamic formulation. Alternatively, one could add pro- 
cess noise to the satellite state itself (position or velocity). 
While this works well in a kinematic solution, it renders the 
satellite states discontinuous between batches. This com- 
plicates incorporation of a deweighted dynamic model in a 
reduced-dynamic formulation and is ill-suited to many sci- 
ence applications in which continuous orbits are required. 

The dynamic, kinematic, and reduced-dynamic tech- 
niques are compared qualitatively in Fig. 1. The dynamic 
solution (top) adjusts the fewest parameters, preserving 
maximum data strength and yielding the lowest formal 
error (error due to data noise), but can suffer a large sys- 
tematic error from mismodeled dynamics. The kinematic 
solution (center) eliminates modeling error, but the or- 
bit transition is determined entirely from the observations, 
data strength is depleted, and the formal error can grow 
large. The reduced-dynamic solution (bottom) optimally 
combines the two techniques to achieve the lowest overall 
error. 

Here, a Kalman filter formulation of the reduced- 
dynamic technique is first presented, then a covariance 
analysis carried out to evaluate its performance in com- 
parison with the dynamic and kinematic approaches is de- 
scribed. Guidelines for selecting a near-optimal reduced- 
dynamic weighting are given, sensitivity to this weighting 
is assessed, and practical aspects of the technique are dis- 
cussed. 


II. Formulation 

The reduced-dynamic technique can be described math- 
ematically in terms of a Kalman sequential filter formula- 
tion. This involves twx> steps: a time update , which makes 
use of a state-transition model to propagate the satellite 
state estimate and covariance from one time batch to the 
next, and a measurement update , which incorporates a new 
batch of measurements. These two steps alternate until all 
data batches are incorporated. 
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A. Time Update 

Let xj be the user satellite state estimate at time tj , us- 
ing data up to the time tj , and x J + i the predicted state es- 
timate at time tj+\ , using data only up to tj \ let j) 

denote the state transition from tj to tj+i- Process noise 
parameters p representing a fictitious 3-D force on the 
user satellite are now introduced. This gives the following 
dynamic state-transition model for the augmented state 
X=[x,p] T and its associated covariance P [8]: 


Xj+ 1 — <f>j Xj + Bwj (1) 

and 

P J+l =4> } Pj4>J + BQ 3 D t ( 2 ) 

where 
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The transition matrix relating Xj+i to the process noise 
parameters Pj is <j> xp (j + l,j); Mj is a 3 x 3 diagonal 
matrix with its ith element 


mi = exp [-(tj+i - tj)/ri] (5) 


B. Measurement Update 

The model for a measurement update in the reduced- 
dynamic technique is the same as in the dynamic or kine- 
matic techniques, with the exception that X and P are 
now associated with the augmented state. Thus, 

X j = Xj+Gj(z j -A j X j ) ( 6 ) 

and 

Pj = Pj — GjAj Pj (7) 

where zj is the measurement vector at time tj ; Aj is the 
matrix of the corresponding measurement partials with 
X ; -; and Gj is the Kalman gain given by 

Gj = Pi Aj {Aj Pi A] + R,)- 1 ( 8 ) 

with Rj being the error covariance of z j. 


The above models have been formulated in terms of cur- 
rent state for clarity. A pseudoepoch state, U-D factorized 
formulation [8] 1 of these models has been implemented in 
the GPS analysis software system, OASIS, developed at 
JPL [9]. 1 

III. Covariance Analysis 

A. Assumptions and Approach 


w j is a white-noise process of covariance Qj> which, for 
convenience, is assumed diagonal with its ith element 
q i — (1 — m; 2 )cTj 2 ; I p is a unit matrix; is the correlation- 
time constant, which controls the decay rate of the corre- 
lation between time batches; and cq is the steady-state un- 
certainty, which is equal to the root-mean-squared (rms) 
value of the process-noise uncertainty after a long time. 
Both <Ti and r, can be selected to be the same for all i in 
this application, so the subscript i will be dropped from 
now on. The relative weighting of the dynamics is var- 
ied by selecting different values for the steady-state un- 
certainty <x, the correlation time r for the process-noise 
parameters, and the a priori uncertainty c r 0 , which is the 
initial error of the parameters. Increasing r and decreasing 
(j Q and a increases the weight on the dynamic information. 
When r — ► oo, a — ■ + 0 and <r 0 — ► 0, the technique reduces to 
conventional dynamic tracking; when r — ► 0,<r — ► oo and 
c t 0 — ► oo, it becomes purely kinematic. It follows that an 
optimal reduced-dynamic solution must be as good as or 
better than both the purely dynamic and purely kinematic 
solutions. 


A covariance analysis comparing the accuracy of 
TOPEX altitude determination expected with the 
reduced-dynamic, dynamic, and kinematic techniques is 
presented next. Assume a constellation of 18 GPS satel- 
lites placed in six orbit planes. GPS P-code pseudo- 
range and carrier phase data are acquired by a receiver on 
TOPEX and by six globally distributed ground sites. Data 
noise, after a 5-minute integration and dual-frequency cor- 
rection for ionospheric delay, is put at 5 cm and 0.5 cm 
for pseudorange and carrier phase, respectively, which is 
consistent with the performance of modern GPS receivers. 
Carrier phase biases are adjusted with a large a priori un- 
certainty. A 2-hour data arc covering a full TOPEX orbit 
is used initially. The ground track of the TOPEX orbit and 
the positions of the six ground sites are shown in Fig. 2. 
Other error sources evaluated are given in Table 1. 


1 S. C. Wu, W. I. Bertiger, J. S. Border, S. M. Lichten, R. F. Sun- 
seri, B. G. Williams, P. J. Wolff, and J. T. Wu, OASIS Mathe- 
matical Description, V. 1.0, JPL D-3139 (internal document), Jet 
Propulsion Laboratory, Pasadena, California, April 1, 1986. 
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The clocks on all GPS satellites and at all but one 
ground site are modeled as white process noise and thus 
are adjusted independently at each time point. This clock 
model is comparable to using doubly difFerenced data and 
is the most conservative (pessimistic) model we can use, 
since it maximizes the number of adjusted parameters 
and hence the formal error. The gravity-error model is 
derived by scaling the difference between corresponding 
coefficients, up to degree and order 20, of two different 
Earth models, GEM- 10 and GEM-L2 [10,11]. This form 
of gravity-error model is convenient to implement, is eas- 
ily varied with a single scaling factor, and has proven re- 
liable in numerous studies over the years. Alternatively, 
one can employ the covariance matrix associated with a 
single Earth model. Our own experiments have shown 
that a 50 percent scaling of GEM- 10 - GEM-L2 approxi- 
mates the error obtained with the covariance matrix from 
the gravity model, GEM-T1 [12], which is one of the best 
current models. 

The 1-cm zenith troposphere error assumes the use of 
a water-vapor radiometer at each ground site. An earlier 
analysis showed that the error from mismodeling atmo- 
spheric drag is less than 1 mm for TOPEX over several- 
hour arcs of data [3]; therefore, drag is not included here. 
Such potential error sources as thermal imbalances and 
outgassing; instrumental delay variations (in the GPS re- 
ceivers) not common to all signals; and imperfect knowl- 
edge of the satellite center of mass, the GPS antenna phase 
center, and platform attitude are being carefully controlled 
for TOPEX and are expected to be below 1 cm. One other 
potentially serious error source is GPS signal multipath. 
Because multipath is not readily treated by covariance 
analysis, separate simulation studies, incorporating all ma- 
jor reflecting surfaces, antenna gain patterns, and receiver 
tracking characteristics have been carried out. The result 
is that while the instantaneous multipath can at times look 
alarming (tens of centimeters on pseudorange and up to 1 
cm on the carrier), it oscillates with periods of minutes or 
less and therefore averages down quickly; the net orbit er- 
ror is typically at the 1-cm level or below after a few hours 
of averaging. 

The inherent GPS data strength allows accurate simul- 
taneous adjustment of GPS orbits and all but three ground 
sites [2-4]. Because the GPS satellites are at an altitude 
where dynamic mismodeling is very small, their solutions 
remain dynamic in all three approaches. Since we know 
that GPS orbits will be routinely tracked with high ac- 
curacy by a global network, a tight a priori uncertainty 
is applied to GPS states. By contrast, to examine the 
strength of each solution technique, a large a priori error 


is placed on TOPEX. In actual operation, a much tighter 
TOPEX a priori error could be used. 

JPL’s recently developed Deep Space Network (DSN) 
Rogue GPS receiver is currently planned for use at all 
ground reference sites. The Rogue can track up to 
eight GPS satellites simultaneously, while the operational 
TOPEX receiver will have a six-satellite capacity. For 
study purposes it is first assumed that the TOPEX and 
ground receivers can observe all visible GPS satellites (typ- 
ically six or seven); then cases with lesser TOPEX receiver 
capacity are examined. The three tracking techniques will 
be assessed here by comparing the TOPEX altitude errors 
over the entire data span. For this, the pseudoepoch-state 
covariances of TOPEX are first smoothed backward and 
then mapped to all time points when data are taken. Com- 
parison is made between the rms errors calculated over the 
entire data span. 

B. Results 

In a preliminary study not shown here, we examined the 
limiting cases of the reduced-dynamic technique. With r 
set to 0 and both <r Q and a set to a large number, the error 
estimate indeed approached the kinematic solution. With 
r set to a large number and both cr 0 and a set to 0. it 
yielded the dynamic estimate. Here, a series of intermedi- 
ate values for r, <r 0 , and a is studied. In general, when r 
is long compared to the batch size, the results vary with 
the batch-to-batch uncertainty <7bb = (1 — m 2 ) 1 / 2 ^, rather 
than with the steady-state uncertainty cr and r individ- 
ually. (The batch-to-batch uncertainty is the one-sigma 
change in value from one time batch to the next that is 
allowed for the process-noise parameters.) In the rest of 
this analysis, a batch size of 5 minutes and a constant 
t — 15 minutes are used for all cases with the reduced- 
dynamic approach; only <t 0 — a is varied to yield a nearly 
optimal solution. 

Figure 3 shows the TOPEX altitude error as a function 
of the percentage of the GEM- 10 — GEM-L2 error, for 
various values of a. Also included are the results with 
dynamic tracking (<r = 0,r — ► oo) and kinematic tracking 
(a — ► oo,r = 0). It is clear that for any finite dynamic- 
model error (in this case dominated by gravity), a range 
of a exists over which TOPEX altitude error is lower than 
with either the dynamic or kinematic solutions. In other 
words, the reduced-dynamic technique is superior provided 
that the dynamic model is properly weighted. A procedure 
for estimating the proper weight is outlined later. 

Figure 4 compares the reduced-dynamic solution with 
the dynamic and kinematic solutions for three different 
viewing capacities for the TOPEX receiver: four, five, and 
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all GPS satellites (typically six, seldom more than seven) 
visible above a 90-deg zenith angle. In the cases with re- 
stricted receiver capacity, the GPS satellites observed are 
carefully selected to minimize satellite switches over the 
observing period (thereby maximizing continuity in carrier 
phase measurements) while still maintaining good observ- 
ing geometry (low PDOP). The gravity error is fixed at 
50 percent of the difference between GEM-10 and 
GEM-L2. A near-optimal weight (a = 0.5 /mi/sec 2 ) is 
used for the reduced-dynamic solution in all three cases. 

When the TOPEX receiver can track all visible GPS 
satellites, geometry is always strong and kinematic track- 
ing is effective; incorporating additional dynamic infor- 
mation through the reduced-dynamic technique improves 
accuracy by only 1 cm. A lower gravity error, perhaps 
achieved through gravity tuning, would of course improve 
reduced-dynamic performance. At the other extreme, 
when the receiver can track only four GPS satellites, ge- 
ometry is often poor and dynamic tracking is far superior 
to kinematic; the optimal reduced-dynamic combination 
again offers little advantage. If, however, the gravity error 
is doubled, as would be the case with a lower orbit, the er- 
ror with dynamic tracking nearly doubles, to 24 cm, while 
the reduced-dynamic performance degrades only moder- 
ately, to 16 cm, illustrating that even weak geometry can 
be of value when dynamics are poorly known. In the inter- 
mediate case, with TOPEX tracking up to five satellites at 
once, the dynamic and kinematic solutions are better bal- 
anced, achieving 12 and 16 cm, respectively. The reduced- 
dynamic combination improves this to 9 cm. In general, 
the reduced-dynamic technique is of greatest value when 
the kinematic and dynamic solutions are comparable. 

Dynamic tracking performance degrades over regions 
where gravity is poorly known (for example, over the 
oceans). Kinematic performance, on the other hand, can 
vary dramatically with changing observing geometry. In 
the reduced-dynamic solution, the two techniques comple- 
ment one another, and the solution is better balanced. 
This is illustrated in Fig. 5, which compares TOPEX al- 
titude accuracy over the whole orbit (two hours), using 
the three techniques. In this case both a TOPEX receiver 
tracking five GPS satellites and the 50 percent of GEM- 
10 — GEM-L2 gravity error are assumed. Both the dy- 
namic and kinematic solutions show peak errors of 25 cm 
or higher at some points. The reduced-dynamic solution, 
with a near-optimal weight ( a = 0.5 /im/sec 2 ), smooths 
these peaks and remains below 13 cm for the entire pe- 
riod. Reduction of the error peaks results from a balance of 
state-transition information between dynamics and kine- 
matics. When the information from one source is weak, 
the Kalman filter places greater weight on the other to 


minimize the overall error. To further illustrate this trade, 
Fig. 6 breaks down the TOPEX altitude error into its con- 
tributing components at the three times (20, 60, and 90 
minutes past epoch) when either the dynamic or the kine- 
matic error grows large. The balance of transition infor- 
mation in the reduced-dynamic solution has resulted in a 
more uniform contribution from all error components. 

In the examples up to this point, a 2-hour tracking arc 
has been used. In general, as the span is increased, the ef- 
fects of data noise and troposphere are reduced while the 
dynamic modeling error grows. In a purely dynamic solu- 
tion, the effect of increasing model error soon overtakes the 
decreasing data error and the overall error tends to grow 
with data span. In the optimal reduced-dynamic solution, 
however, the estimator continuously shifts weight to the 
increasingly strong data, away from dynamics, as the span 
increases. This deweighting of dynamics is a natural conse- 
quence of the estimation process; no change in a is needed 
since the optimal cr applies to a specific dynamic-model 
error, independent of data span. As a result, with optimal 
weighting, the overall performance will tend to improve 
with increased data span. 

Figure 7 compares the TOPEX altitude error when 
2-hour and 4-hour data spans are used. The longer data 
span reduces the error over the initial 2-hour period to 
7 cm from the 8.9 cm obtained with the original 2-hour 
arc. An examination of the error breakdown shows a re- 
duction in gravity error, reflecting the deweighted dynam- 
ics, as well as in other errors. Although spans longer than 
four hours have not been studied for TOPEX, we expect 
that the error will reduce monotonically with data span. 
Because the weight on the dynamic model decreases with 
longer data span, a reduced-dynamic solution will tend to 
a kinematic solution as the span is increased. Note that 
this is true only on the assumption that a fixed dynamic 
model is used, independent of the data arc length. If the 
model is improved through tuning or other efforts, the op- 
timal weight for a given data span will shift back toward 
the dynamic solution [13]. 

C. Other Applications 

TOPEX nicely illustrates the benefits of reduced- 
dynamic tracking since its altitude of more than 1,300 km, 
six-satellite receiver capacity, and relatively compact di- 
mensions permit both good observing geometry and rea- 
sonably well-modeled dynamics. A far greater modeling 
challenge is presented by several other current or planned 
NASA space platforms: the large (14 m) platforms of the 
polar-orbiting Earth Observing System (EOS), which will 
carry heavy slewing instruments and fly at 700 km; the 
actively maneuvering space shuttle, at altitudes as low as 
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300 km; and the sprawling (155 m) Space Station Freedom, 
at about 400 km. All will eventually carry experiments 
seeking tracking accuracies better than 10 cm. Indeed, a 
recent international workshop on space geodesy set a goal 
of “no more than 1 cm rms error, single pass, without or- 
bit discontinuities” [14] for tracking future orbiting ocean 
altimeters, such as the one that will fly on EOS. 

Since one cannot expect to approach centimeter or even 
decimeter accuracy in modeling the dynamics of such un- 
gainly platforms, the optimal orbit solution strategy will 
be almost purely kinematic. To maximize performance un- 
der kinematic tracking, geometric observing strength must 
be maximized. With that in mind, one can take the ex- 
amples of EOS (98-deg inclination, 705-km altitude) and 
Space Station Freedom (28 deg, 400 km) and carry out co- 
variance studies under a more robust set of assumptions: 
the GPS constellation is increased to 24 satellites, as is ex- 
pected to occur by 1995; the ground network is expanded 
to 10 sites; the flight receiver is extended to track all satel- 
lites in view down to the Earth limb (typically a dozen 
or more); and the three fixed ground sites are assumed 
known to 2 cm in each component, which is expected to 
be achieved or surpassed by very long baseline interferom- 
etry within the next few years [14]. 

Figure 8 plots the rms position errors estimated for a 
purely kinematic solution for arc lengths of 2, 4, and 8 hrs. 
Under the revised assumptions, the observing geometry 
is so consistently strong that few-centimeter accuracy is 
achieved continuously, and the rms error approaches 1 cm 
per component after 8 hrs. Despite the widely different or- 
bits and dynamics, performance is virtually indistinguish- 
able between EOS and the Space Station, reflecting the 
full emphasis on geometric strength (which differs negligi- 
bly between the two) over dynamics. Note that the purely 
kinematic position solution is referred to the phase center 
of the orbiter’s GPS antenna, which can be calibrated with 
submillimeter accuracy, rather than to the platform center 
of mass, which can be difficult to pinpoint on a large and 
variable structure. Because the dynamics of these orbiters 
are so poorly known, optimal reduced-dynamic solutions 
would improve the kinematic results by only 1 or 2 mm — 
at a great cost in computation. 

For satellites in high-Earth elliptical orbits, which 
might cover a range of altitudes from 1,000 km to 
20,000 km or higher, a variety of different process-noise 
force models for reduced-dynamic tracking can be em- 
ployed. At the lower altitudes, models similar to those 
adopted for TOPEX can be utilized, while at the high alti- 
tudes (where gravity and drag are insignificant), low levels 
of process noise with long time constants are more appro- 
priate to absorb unmodeled accelerations that are due to 


gas leaks and solar radiation pressure. Covariance analyses 
to study orbit determination performance for high-Earth 
elliptical orbiters using GPS will be reported in a future 
article; however, preliminary analyses indicate that orbit 
accuracies at the several-decimeter level should be possible 
even at altitudes of 20,000 km or higher. 2 

IV. Weighting the Dynamic Model 

For applications such as TOPEX where geometry and 
dynamics are more balanced, a procedure is needed to es- 
timate the weight for the dynamic model, specified by 
cr with any adopted r, that minimizes the rms orbit er- 
ror. This may be difficult to do precisely, since the quality 
of dynamic models is not always well understood; often, 
in fact, unsuspected modeling errors are present. Fortu- 
nately, the sensitivity to a departure from optimal weight- 
ing appears to be low. This is illustrated in Fig. 9, 
in which the TOPEX altitude error with the reduced- 
dynamic technique is shown as a function of the level of 
gravity error for three different weightings. Dynamic and 
kinematic results are shown for comparison. The weight 
cr — 0.5 /tm/sec 2 is nearly optimal for a gravity error of 
50 percent of GEM- 10 — GEM-L2. The two other curves 
are for values of a a factor of 4 larger (<7 = 2 /jm/sec 2 ) and 
smaller (<r = 0.125 /i m/sec 2 ). This wide range of subopti- 
mal weights increases the TOPEX altitude error by only 
0.3 cm at one end and 1.5 cm at the other. In other words, 
the performance is fairly insensitive to suboptimal a. 

In practice, an approximate weight can be estimated in 
advance through a covariance analysis. This is done by 
a series of filter runs, each using a different weight, sim- 
ulating the actual measurement and estimation scenario. 
Realistic data noise and models for all significant (includ- 
ing dynamic) error sources must be considered and their 
effects on the orbit determination evaluated. The weight 
resulting in the lowest estimated rms orbit error is the best 
estimate of the optimal weight for actual data processing. 

A misjudgment of the dynamic-model error will, of 
course, result in selection of a suboptimal weight; how- 
ever, one can take care to minimize the effect of such a 
misjudgment with the following strategy: The assumed 
dynamic error model is used first to predict the perfor- 
mance of both the dynamic and kinematic solutions. If 
either of these appears far superior to the other, say, by 
a factor of 3 or more, the slight improvement that would 


2 S. M. Lichten, “Orbit Determination for High-Earth Elliptical Or- 
biters Using GPS I. Initial Results for VSOP,” JPL Section 335 
IOM 335.4-89-123 (internal document), Jet Propulsion Laboratory, 
Pasadena, California, November 15, 1989. 
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result by combining the two approaches with the reduced- 
dynamic technique may not justify the extra effort, and the 
simpler form can be adopted. If neither technique is far 
superior, a weight departing from the predicted optimum 
in a direction favoring the kinematic (i.e., larger <r) will be 
prudent. Such a bias in the weighting can reduce the more 
damaging effect of dynamic error if the error is larger than 
expected. This point is illustrated in Table 2, which sum- 
marizes the effect of weight misjudgment by a factor of 4 
on TOPEX altitude determination for two different levels 
of gravity error. These results suggest that a weight biased 
in favor of the kinematic approach is preferable when the 
level of dynamic error cannot be well determined. 

V. Other Considerations 

In the analysis of the reduced-dynamic solution, a fic- 
titious 3-D force on TOPEX was treated as process noise 
and adjusted together with TOPEX and GPS states. In- 
troduction of the force is merely for the purpose of chang- 
ing the filter model to reduce that modePs reliance on the 
dynamic model. Since, in the real world, this force does 
not exist, its presence in the formulation adds an error 
source in the estimation process, causing the formal error 
to be overestimated. To remove this effect, an evaluation 
run 3 of the filter is needed. In an evaluation run, the fil- 
ter model would be specified as before, including the ficti- 
tious force, but the contribution of this force in the truth 
model would be ignored. Such an algorithm is fairly com- 
plicated when a smoothing process is required, because of 
the dynamic process-noise parameters involved. For a fair 
estimate of this effect, evaluation runs with the process- 
noise force replaced by correlated piecewise-constant (in 
time) forces have been made, thus avoiding the need of 
smoothing. These runs show that the spurious increase in 
the formal error due to the fictitious force is only a few 
millimeters; the corresponding increase in the total error 
is even smaller, typically 1 to 2 mm. 


VI. Conclusions 

A reduced-dynamic technique for determining the or- 
bits of Earth satellites is made possible with observa- 
tions of the Global Positioning System. In this technique, 
satellite-state-transition information obtained from both a 
dynamic model and continuous GPS carrier phase obser- 
vations is optimally combined to improve orbit determi- 
nation accuracy. Analysis indicates that a significant im- 


3 C. L. Thornton, “Triangular Covariance Factorizations for Kalinan 
Filtering,” TM 33-798 (internal document), Jet Propulsion Labo- 
ratory, Pasadena, California, October 1976. 


provement can be expected when neither of the two types 
of state transition information is far superior to the other. 
Performance is not highly sensitive to the relative weight- 
ing between dynamic and kinematic information. When 
the actual level of dynamic model error is uncertain, an 
additional deweighting of the dynamic model is recom- 
mended; this would prevent an inordinately large error 
resulting from larger-than-expected dynamic model error. 

Although a tracking arc of 2 hours was used for most 
of the reduced-dynamic analysis, a 4-hour span was exam- 
ined to illustrate the improved performance with increased 
data. Further improvement can be expected with longer 
data spans due to reduction of the effects of data noise 
and random error in tropospheric delay modeling. The ef- 
fects of increasing dynamic-model error will be automati- 
cally controlled in the reduced-dynamic solution by further 
deweighting and will thus remain at a low level. Because 
the weight on the dynamic model is lowered with grow- 
ing data strength, a reduced-dynamic solution will grad- 
ually approach the kinematic solution as the data span 
increases, provided that a fixed dynamic model is used. 

Reduced-dynamic tracking can be used with any Earth 
satellite that can adequately observe GPS. The altitude 
range over which reduced-dynamic tracking provides useful 
improvement over dynamic and kinematic tracking will de- 
pend on the actual level of the dynamic-model uncertainty. 
It is expected that satellites at altitudes between 400 and 
2,000 km will receive the greatest benefit. Above this 
range, GPS observability diminishes while dynamic-model 
errors decrease markedly, so greater reliance on dynamic 
tracking will be favored. Below this range, uncertainties 
in gravity and atmospheric drag become so great that a 
kinematic solution may be favored for simplicity. For some 
exceptional vehicles, other considerations apply. The ac- 
tively maneuvering space shuttle may receive no benefit 
from dynamics at any altitude, while drag-compensated 
satellites may exploit dynamics at even the lowest orbit 
altitudes. 

In designing a GPS-based precise tracking system for an 
Earth orbiter, there is a simple trade-off between modeling 
accuracy and geometric strength to consider. Where the 
models are strong, the geometry can be relaxed and the 
flight and ground systems kept relatively simple. Where 
the models are weak, as will be the case with a number of 
dynamically complex missions in the future, the geometry 
must be strengthened. In any case, the global coverage and 
unique mix of data types offered by GPS ensure that there 
will be a practical system design and solution strategy that 
can deliver orbit accuracies well under a decimeter for any 
low-Earth satellite. 
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Table 1. Error model and other assumptions used in 
covariance analysis 


Model component Assumption 


User satellite 

TOPEX (1,334 km in altitude) 

Number of stations 

6 (cf. Fig. 2) 

Number of GPS satellites 

18 

Cut-ofT elevation 

10 deg at stations 
0 deg at TOPEX 

Data type 

P-code pseudorange 
Carrier phase 

Data span 

2 hours 

Data interval 

5 minutes 

Data noise 

5 cm (pseudorange) 
0.5 cm (carrier phase) 

Carrier phase bias 

10 km (adjusted) 

Clock bias 

3 (.isec (adjusted as white 
process noise) 

TOPEX epoch state 

2 km; 2 m/sec (adjusted) 

GPS epoch states 

2 m; 0.2 mm/ sec (adjusted) 

Station location 

5 cm each component 

Zenith troposphere 

1 cm 

Earth’s GM 

1 part in 10 8 

Gravity 

Scaled GEM-10-GEM-L2 
(see text) 

Solar pressure 

10 percent 


Table 2. Effects of weight misjudgment on TOPEX altitude 
determination accuracy 


Gravity error 

x optimal a 

optimal a 

4 X optimal cr 

50 percent 
100 percent 

10.4 cm 
12.3 cm 

8.9 cm 
9.7 cm 

9.2 cm 
10.1 cm 
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TECHNIQUES 


Fig. 3. Performance of reduced-dynamic technique with different Fig. 4. Performance of dynamic, reduced-dynamic, and 

weighting on dynamic model. kinematic techniques. 
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Fig. 5. Variation of TOPEX altitude error over time. 
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Fig. 6. Breakdown of TOPEX altitude error at three different epochs. 
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Fig. 7. Comparison of reduced-dynamic tracking performance with different data spans. 




Fig. 8. Estimated kinematic tracking accuracies for EOS and 
Space Station, both with a robust GPS observing system. 



Fig. 9. Insensitivity of TOPEX accuracy to suboptimal weighting. 
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